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1 Introduction 



One-sided and two-sided exit problems for the compound Poisson processes and jump 
diffusion processes with two-sided jumps have been apphed widely in a variety of fields. 
For example, in the theory of actuarial mathematics, the problem of first exit from a half- 
line is of fundamental interest with regard to the classical ruin problem and the expected 
discounted penalty function or the Gerber-Shiu function as well as the expected total 
discounted dividends up to ruin. See e.g. Kliippelberg et al. (2004), Mordecki (2004), 
Xing et al. (2008), Cai et al. (2009), Zhang et al. (2010), Chi (2010) and Chi and Lin 
(2011). In the setting of mathematical finance, the first passage time play a crucial role for 
the pricing of many path-dependent options, American-type and Russian-type options, 
see e.g. Kou (2002), Kou and Wang (2003, 2004), Asmussen et al. (2004), Levendorskii 
(2004), Ahh and Kyprianou (2005), Cai et al. (2009), Cai and Kou (2011) and as well 
as certain credit risk models, see, for example, Hilberink and Rogers (2002), Le Courtois 
and Quittard-Pinon (2006) and Dong et al. (2011). Many optimal stopping strategies 
also turn out to boil down the first passage problem for jump diffusion processes, see 
e.g. Mordecki (2002). In queueing theory one-sided and two-sided first-exit problems 
for the compound Poisson processes and jump diffusion processes with two-sided jumps 
have been playing a central role in a single-server queueing system with random workload 
removal, see e.g. Perry et al. (2002). Usually, when we study the first passage problem, 
the models with two-sided jumps are more difficult to handle than those with one-sided 
jumps, because the undershoot and overshoot problem could not be avoided. Despite 
the maturity of this field of study it is surprising to note that, until very recently, it 
can only be solved for certain kinds of jump distributions, such as the Kou's double 
exponential jump diffusion model (see Kou (2002), Kou and Wang (2003)). Recently, Cai 
and Kou (2011) proposed a mixed-exponential jump diffusion process to model the asset 
return and found an expression for the joint distribution of the first passage time and the 
overshoot for a mixed-exponential jump diffusion process. However, determination of the 
coefficients in their expression still remains a mathematical and computational challenge; 
See also Cai (2009) for the case of a hyper-exponential jump diffusion process. In this 
paper we will further study the first passage problems in Cai and Kou (2011) and give an 
explicit expression for the joint distribution of the first passage time and the overshoot 
for a mixed-exponential jump process with or without a diffusion. Moreover, we present 
several applications in insurance risk theory and in finance. 
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The rest of the paper is organized as follows. In Section 2, the model assumptions are 
formulated. In Section 3, we study the one-sided passage problem from below or above 
for compound Poisson process and jump diffusion process. In Section 4, we give explicit 
expression of the Gerber-Shiu function with two-sided jumps. In Section 5, we present the 
analytical solutions to the pricing problem of one barrier options and lookback options, 
and in the last section we derive a closed-form expression for the price of the zero-coupon 
bond. 



2 Mathematical model 



A jump diffusion process X — {X{t) : i > 0} is defined as 

Nt 

X{t)=nt + aWt + Y."^i, (2.1) 

where {Wt, i > 0} is a standard Brownian motion with Wq — 0, {Nt, i > 0} is a Poisson 
process with rate A, constants e M, cr > represent the drift and the volatility of the 
diffusion part respectively, and the jump sizes {Y^; i > 1} are independent and identically 
distributed random variables whose probability density function is given by More- 
over, it is assumed that {Wt}, {Nt} and {Fj} are independent. When a — the process 
(2.1) is the so called the compound Poisson process with positive and negative jumps and 
linear deterministic decrease or increase between jumps according to < or > 0. The 
processes cover many models appearing in the literature such as the compound Poisson 
risk models, the perturbed compound Poisson risk models, and their dual models. From 
now on, we shall denoted by {P^ : x G M} probabilities such that under Pj,, X(0) = x 
with probability one. Moreover, E,j. will be the expectation operator associated to P^- 
For convenience, we shall write P = Pq and E = Eq. 

It is easy to see that X is a special case of Levy processes with two-sided jumps, whose 

infinitesimal generator of X{t) is given by 

1 f°° 
Cg{x) = -a^g"{x) + iig\x) + Xj {g{x + y)- g{x))fY{y)dy, 

for any twice continuously differentiable function g. The moment generating function of 
X{t) is S(e^^(*)) = e'^(^)*, t > Q,U{z) = 0, where '0(z), called the exponent of the Levy 
process X, is defined as 

^(z)^la^z''+i^z + X{E[e'^]-l). (2.2) 
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For more about general Levy processes, we refer to Bertion (1996), Kyprianou (2006) and 
Doney (2007). 

3 First passage problems 

We now turn to one-sided passage problems for the Levy process (2.1). For two flat 
barriers h and H {h < H), dcflne the first downward passage time under h and the first 
upward passage time over H by 

:= inf{t > : X{t) < h}, := inf{t > : X{t) > H}, 

with the convention that inf = oo. In the next two subsections we will investigate 
the distributions of the following quantities: first upward passage time and overshoot 
X{t^) — H; first downward passage time and undershoot h — X(t^). 

3.1 One sided exit from above 

In this subsection we assume that the downward jumps have an arbitrary distribution with 
density /_ and Laplace transform /_, while the upward jumps are mixed-exponential, i.e. 

m 

fyiy) = pf-{-y)l{y<o} + q^Pirjie~'''n{y>oy, (3.1) 

where constants p,q > 0,p + q = 1, < r]i < rj2 <■■■< r]m < oo and Yl^iPi — ^■ 
The Levy exponent of X is given by 

Mz) = \<y'z' + + A f; ^ +pf-{-z) - 1) . 

Using the same argument as in Cai and Kou (2011) we have the following 

Lemma 3.1. (i) For sufficiently large a > 0, if a > or /i > and a — 0, then the 
equation ipi{z) — a has exactly m + 1 distinct positive roots f3i, - ■ ■ , f3m+i satisfying 

< /Si < /S2 < • • • < /3m+i < 00. 
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(a) ///i < and (J = 0, then the equation '4'2{z) = a has exactly m distinct positive roots 
A, • • • ,/3m satisfying 

< Pi < /32 < ■ ■ ■ < Pm < oo. 



Cai and Kou (2011) found the joint distribution of the first passage time and 
X{Tff) in the case a > under the additional assumption f-{y) is also mixed-exponential. 
However, for a general f-{y) as in the case of the upward jumps are mixed-exponential 
(cf. Yin, Shen and Wen (2013)), for any sufficiently large a > 0, ^ < ?7i and x < H, we 
have 



m+l 
k=l 



(3.2) 



where w :— {wi, • • • , Wm+i)' is a vector unique determined by the following system ABw — 
J, here A is an (m + 1) x (m + 1) matrix 



A 



1 



1 

m 



1 



Vl-fil V1-P2 



'Qm Vm 



Vm-fil Vm-/32 



Vl-Pr, 



Vm 
rim-l3m + l 



B is an {m + l) x {m + l) diagonal matrix and J is an (m -|- 1) -dimensional vector 



In this paper we will determined the coefficients w^'s explicitly. Moreover, we also 
consider the cases /i > 0,a = and /i < 0,a = 0. 

Theorem 3.1. For any sufficiently large a > 0, we have 



(i) for 9 < Tji and x < H, 

-art+OX 



( e 



N 



fe=l lii=ll-^ r)i) 



(ii) for y > 0,x < H , 



(3.3) 



(e-"^^l{x^^.Hedy}^ ^ (^^koSo{y) + J2Akne-^^^^ e'^^^^'^^dy, 



(3.4) 
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(Hi) for X < H , 



(3.5) 



(iv) for X < H , y > 0, 



-Ex ( e "'^"1{X +-H>y} 



N 



k=i \i=i 



kie 1 e-^'=(-^-^) 



(3.6) 



(v) for X < H, 



N 



■h{H-x) 



(3.7) 



fe=i 



where Pi, - • ■ , (3^ are the positive roots of the equation — a, 5q{x) is the Dirac delta 

at X — and 

m + 1, if (7 > 0, or (7 = and > 0, 
m, if cr = and fJ> < 0, 



N 



mii(i-S) 
nfc=i,fc^j(i - ^) 



fcO 



^JM^, if (7 > 0, or (7 = and > 0, 



0, 



if (7 = and < 0, 



-^kl — T-rm 7 ' — ■■ ■ ) 



Proof We prove the result for the case a > only, the rest cases can be proved 
similarly. To prove Theorem 3.1, the most difficult part is to find the inverse of matrix 
A. For simplicity, we write 

' All Ai2 
A21 A22 



A^ 



where 



and 



Al = (l), A12 = (!,••• ,l)ixm, ^1 = 



m-P2 



Vl 



A 



Vn 



Vl- Pi' ' Vm- Pi 



22 



L ■nm-P2 



rim—Pm+l 
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Note that A22 can be written as A22 = JiCi, where Ji = Diag{rii, ■ ■ ■ ,1]^} is a 
diagonal matrix, Ci = { ^._^.^^ }i<i,j<m is a Cauchy matrix of order m which is invertible 
and the inverse is given by Cf ^ = {dij}mxm, where 



dij = {r]j-f3i+i) — 



Here 



i=l 



Then the inverse of A22 is given by 



^22 — 



m 



21 



2m 



ml 



The determinant of Ci is given by (see Calvetti and Reichel (1996)): 

Ul<i<j<rrXVi - Vj)il3j+l - A+l) 



det(Ci) = 



After some algebra, 



where 



A/A 



22 



nr=i(A+i-/9i) 



A/A22 := An - Ai2A2^A2i 



is the Schur complement of the block A22 in A, which is a matrix of order 1. By Schur's 
formula (see Zhang (2005)), 

det(A) = det(A22) ■ det{A/A22) 0. 

Moreover, by Banachiewicz inversion formula (see Zhang (2005)), the inverse of A is given 

by 

= {A/A22)-' \ \ , , 

[ ->l22^21 ^22 ^21^12^22 +^22 (^/^22) 

After some algebra, we have 



j-M A 4-1-1- 4-11-4/4 ^-^ ^ife) V I nr.i(ft+i-ft) ^ 

A,, A,,A,.A^ + A^ (A/ A,,) - j^-^ ^ — + -^^^^—^^d,^ 



l<«j<m 



Now by solving ABw — J we find that 



/ 77'"+^ (1 - ^) Tjm+l (I - ±] \ 

r ni.(i-i) • • ni>(i-i) )• 

from wliicli and (3.2) we get (3.3). 
By tfie fractional expansion, 

— -— — ^-Afeo + ^W- ^^fcm- ^, (3.8j 

where the coefficients Afej's are defined in the theorem. Substituting (3.8) into (3.3) and 
inverting it on 9 immediately lead to (3.4). (3.5)-(3.7) are direct consequence of (3.4). 
This ends the proof of Theorem 3.1. 

Example 2.1 Let m — 1, several expressions obtained by Theorem 3.1. When a > 
or a — and /j, > 0, ior x < H,9 < rji and y > 0, we recover the following three formulae 
which obtained by Kou and Wang (2003): 

When a — and // < 0, then ior x < H,9 < rji and y >0, 



V / r]i-e 

E (p-^^H^ + \ - p-my^]lIL^p-0i{H-x) 



3.2 One sided exit from below 

In this subsection we assume that the upward jumps have an arbitrary distribution with 
Laplace transform /_|_, while the downward jumps are mixed-exponential, i.e. 

m 

fY{y) = phiy) + qJ2p,Vje'>^n{y^o}, (3.9) 



where constants p,q>0,p + q = l,0<rii<ri2<---<rim<oo and jyjLiPj — 1- -^y 
(2.2), the Levy exponent of X is given by 



By replacing X by —X in the previous section, we get the main finding in this section. 

Theorem 3.2. For any sufficiently large a > 0, we have 
(i) for 9 > Q,x > h, 



llt=l,t^feV rj _rk{x-h) (3 10) 



(ii) for X > h,y > 0, 

£;(^e-"^'ri{,_x^_erf,}^ -JlB, (^A,oSo{y) + J2Ame-'^'''^ e-^'-^^-^^dy, (3.11) 
(in) for X > h, 

(e-"^^l{x .=h}) = Yl BkAkoe-'^'^''-'\ (3.12) 

(iv) for X > h, 

(v) for X > h, 

J 

{e-^^^) = Bke-'^'^''-^\ (3.14) 
fc=i 

where — ri, • • • , — rj are the negative roots of the equation "02 (^) = Oi, and 

J m + 1, (7>0, or(7 = and < 0, 

m, (J = and > 0, 



)J / m \ J 

= E E ^'^O '^"''^''"'^ = E ^'^ (1 - ^^o) e-^'=(^-'^\ (3.13) 
fc=l \z=i / k=l 



J'^^''' , (T > 0, or (7 = and // > 0, 
0, 0" and < 0, 



9 



Remark 3.1. The result (3.14) o,grees with the result of Theorem 1.1 in Mordecki (2004), 
where only the case a > and Pi > {i = 1, ■ ■ ■ ,m) is considered. 

Example 2.2 Letting m = 1 in Theorem 3.2. When cr > or cr = and < 0, for 

9 < Tji and y >0, 

—r2{x—h) \ 

V(^ + ^i)(^2-ri)' ^ {9 + r|,){r,-r^f J' 

E (p-'^^h-i \ - ^-m/ (^2 -?7l)(?7l ( -ri{x-h) _ -r2{x-h)\ 

EJe-^'^h) = ^2fa ^ ri(r2 -?7i) ^_^^(^_fc)^ 

?7i(r2-ri) rji(f2-ri) 
When (7 = and A* > 0, then for 9 < r]i and y > 0, 

V / 9 + r]i 

F. fp~°''^h-l \ — p-Viy'!]l^l^p-riix-h) 

^x\e -e e 



4 Applications to Gerber-Shiu functions 

We consider an insurance risk model in which the insurer's surplus process is defined as 

Nt 

U{t) =u + X{t) ■.= u + nt + aWt + Y^ Yi, t > 0, (4.1) 

i=l 

where X{t) is defined by (2.1) with jump density (3.9). The time of (ultimate) ruin is 
defined as r = inf{t > : Ut < 0}, where r = oo if ruin does not occur in finite time. As 
applications, we obtain the following special case of the Gerber-Shiu functions for surplus 
processes with two-sided jumps. 

0(ii) = E(e-''^w(\U(T)\)l(T < oo)|C/(0) = u), 

M^) = E{e-''^w{\U{r)\)l{r < oo, U{t) = 0)|C/(0) = u), 

0,(«) = E{e-''-w{\U{T)\)l{T < oo, U{t) < 0)|t/(0) = u), 

where a > is interpreted as the force of interest and w is a non-negative function 
defined on (— oo, 0]. Note that a more general form of Gerber-Shiu function was originally 
introduced in Gerber and Shiu (1998) for the classical risk model. 

Prom Theorem 3.2 (ii) we get the following result. 
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Corollary 4.1. Suppose that U{t) drifts to +00, then we have 

POO 

<j>{u) = / w{y)Ki^\y)dy, (4.2) 
J 

Uu)^w{0)J2BkAkoe-'''\ (4.3) 

k=l 

J / m „oo \ 

Uu) = E E / Hy)e-'"'dy e-^'\ (4.4) 
k=i \i=i -^0 / 

where 's, A^i 's and 's are defined as in Theorem 3.2, and 

J / m \ 



fe=l \ 1=1 



Remcirk 4.1. We compare our results with the existing literature. For the case of a — 
andY has a double exponential distribution, the result (4-2) was found by Cai et al (2009); 
For (7 = and fj, — 0, the result (4-2) was found by Albrecher et al. (2011, see (3.2)); 
For II — 0, the result (4-2) was found by Albrecher et al. (2011, see (9.3)); For (7 = 
and n < 0, the results (4.2)-(4.4) were found by Eric C. K. Cheung (see Albrecher et al. 

(2011, PP.443-444))- 



5 Applications to pricing path-dependent options 



As applications of our model in finance, we will study the risk-neutral price of barrier and 
lookback options. These options have a fixed maturity T and a payoff that depends on the 
maximum (or minimum) of the asset price on [0, T]. The asset price process {S{t) : t > 0} 
under a risk-neutral probability measure P is assumed as S{t) = S'(0)e^'^*\ where X{t) is 
given by (2.1), S'(O) = e^^^\ We are going to derive pricing formulae for standard single 
barrier options and lookback options, based on the results obtained in Section 3. 



5.1 Lookback options 

The value of a lookback option depends on the maximum or minimum of the stock price 
over the entire hfe span of the option. Let the risk-free interest rate be r > 0. Given a 
strike price K and the maturity T, it is well-known that (see e.g. Schoutens (2006)) using 
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risk-neutral valuation and after choosing an equivalent martingale measure P the initial 
(i.e. t = 0) price of a fixed-strike lookback put option is given by 

Lj;,(i^, T) = e-'-^E ( sup S{t) - kY ; 

\0<t<T ) 

The initial price of a fixed-strike lookback call option is given by 

The initial price of a floating-strike lookback put option is given by 

^/;oa*.n.(r) = ( sup S{t) -sX ■ 

\0<t<T J 



The initial price of a floating-strike lookback call option is given by 

-c / 

-'floating \ 



L^aUn,iT) = e-^^E(ST-MS{t) 



0<t<T 

In the standard Black-Scholes setting, closed-form solutions for lookback options have 
been derived by Merton (1973) and Goldman et al. (1979). For the double mixed- 
exponential jump diffusion model, Cai and Kou (2011) derived the Laplace transforms of 
the lookback put option price with respect to the maturity T, however, the coefficients 
do not determinate explicitly. 

We shall only consider lookback put options because lookback call options can be 
obtained similarly. For jump diffusion process (2.1) with jump size density (3.1), the 
condition 771 > 1 is imposed to ensure that the expectation of e~'^*St well defined. 

Theorem 5.1. For all sufficiently large 5 > 
(i)JorK>So, 

nf -I Pi,r + S \ 



[II 

i=l 

where (3i^r+s, ' ' ' , I^N,r+s o,re the N positive roots of the equation ipi{z) — r + 5, and 



/ -STtP (rr\jrp _ ^0 1 = 1 ^ ^ I ^0 ^0 



N 



m + 1, (7 > 0, or (7 = and A* > 0, 
m, (7 = and < 0. 
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Proof (i). We prove it along the same line as in Cai and Kou (2011). Set A; = In ^ > 0, 
then 



-rT 



e^P ( sup X{s)>y] dy. 

yO<s<T 



It follows that 

re-^^Lf^,{K,T)dT = Sorey[re-^^^'^^F{sup,,,,rX{s)>y)dT]dy 



The result follows from Theorem 3.1 and (5.1). 
(ii). Since 



exp I sup X(t) 

\0<t<T 



— So, 



it follows that 

r e-'^Lj,^^un,{T)dT = So J.^ e-(^+^)^E [exp (supo<,<^ X{t))] dT - f 

= ^^jE [exp (supo<,<e(,+5) Xit))] - f 

1 + ir (suPo<s<e(r+5) X{s) > y) dy] - 
"l + /~e^E(e-('-+^)^-^)" 

The result follows from Theorem 3.1 and (5.2). 



Sa 
r+S 

So 
r+S 



So 



So 

s ■ 



(5.1) 



(5.2) 



5.2 Barrier options 

The generic term barrier options refers to the class of options whose payoff depends on 
whether or not the underlying prices hit a prespecified barrier during the options' lifetimes. 
There are eight types of (one dimensional, single) barrier options: up (down)-and-in (out) 
call (put) options. For more details, we refer the reader to Schoutens (2006). Kou and 
Wang (2004) obtain closed-form price of up-and-in call barrier option under a double 
exponential jump diffusion model; Cai and Kou (2011) obtain closed-form expressions of 
the up-and-in call barrier option under a double mixed-exponential jump diffusion model. 
Here, we only illustrate how to deal with the down-and-out call barrier option because the 
other seven barrier options can be priced similarly. For jump diffusion process (2.1) with 
jump size density (3.9), given a strike price K and a barrier level U, under the risk-neutral 
probability measure P, the price of down-and-out call option is defined as 

DOC = exp(-rr)E[(5T - i^)+l(i„fo<,<,5*>f/)|5'o], U < S^. 
Let /i = In I- and A; = - In ir. Then 

So 

DOC{k,T) := DOC = exp(-rr)E^[(5oe^^ - e-'')+'^(r->T)]- 
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Theorem 5.2. For any < 4> < r]i — 1 and r + > + 1), then 



r-oo poo 



e-'t>k-^'^DOC{k,T)dkdT 



where where —Ri, ■ ■ ■ , —Rj are the negative roots of the equation ip2{T) = r + tp, and 



J 



m + 1, (7 > 0, or (7 = and < 0, 
m, a — and > 0, 



Proof Using the same argument as that of the proof of Theorem 5.2 in Cai and Kou 
(2011), we get 



joo joo^ e-'^k-^^DOC{K T)dkdT = J^^ JZ, e-^'^-('-+^)^E.[(5oe^- - e-'^yi^,-^^^]dkdT 



0(</>+l) <p+r 

and the result follows from Theorem 3.2 (i). 



6 The price of the zero-coupon bond 

In this section, we give a simple application on the price of the zero-coupon bond under 
a structural credit risk model with jumps. As in Dong et al. (2011), we assume that the 
total market value of a firm under the pricing probability measure P is given by 

V{t)^Voe^^*\ t > 0, 

where Vq is positive constant, X{t) is defined as (2.1). For > 0, set m = — ln(i^/Vo), 
define the default time as 

Tu = inf{i : V(t) <K} = mi{t :u + X{t)<0}. 

Given T > and a short constant rate of interest r > 0, Dong et al. (2011) shown that 
the Laplace transform of the fair price B{0, T) of a defaultable zero-coupon bound at time 
with maturity T is given by 

^, , 1 - £;[e-(^+-)-"] i?£;[e-(^+^)-«F(r„)l(r„ < oo)] 

= \ \ 7? > 

7 -I- r A7 
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where i? G [0, 1] is a constant. When the jump size distribution is a double hyperex- 
ponential distribution, a closed-form expression is obtained, but the coefficients can not 
determined explicitly (except for n = 2). Now applying the result in Section 3.2, we get 
the following result: 

Corollciry 6.1. If the process X{t) is defined as (2.1) has jump size density (3.9), we 
have 

7+^ 7^ ni=i(i + i) 

where —pi, • • • , —pj are the negative roots of the equation ■02 (p) = 7 + and 



J 



m + 1, (J > 0, or (J = and A* < 0, 

m, 0" = and /i > 0, 

T-rm Q _ PJ_^ 

C,= \^'='^ , J = !,■■■, J. 
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